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Abstract 

Entangled states of two ions are realized by using an adiabatic process. Based on the proposal by Lin- 
ington and Vitanov | ij], we have generated Dicke states in optical qubits of two ^^Ca.^ ions by applying 
frequency-chirped optical pulses with time-dependent envelopes to perform rapid adiabatic passage on side- 
band transitions. One of the biggest advantages of adiabatic approaches is their robustness against variations 
in experimental parameters, which is verified by performing experiments for different pulse widths or peak 
Rabi frequencies. Fidelities exceeding 0.5, which is the threshold for inseparable states, are obtained over 
wide ranges of parameter values. 

PACS numbers: 03.67.Bg, 37.10.Ty 
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I. INTRODUCTION 



Multipartite entanglement plays a central role in quantum information science. Entangled states 
have been generated and studied in various systems, including photons and trapped ions [l2|-[6|]. Of 
multipartite entangled states, a set of states known as Dicke states [|7D has recently been attracting 
interest. Dicke states are entangled states that are symmetric with respect to particle permutations. 
They are defined as 



tn 



-1/2 



(1) 



where A'^ and m are respectively the number of particles and excitations and YikPki---) indicates 
the sum over all particle permutations. W states, a subset of Dicke states with a single excitation, 
have been produced and investigated using up to eight trapped ions [j2|]. Dicke states of up to six 
photons have been generated and their unique characteristics, such as entanglement persistency 
against photon loss and projective measurements, have been studied DSilalsilSD. 

Linington and Vitanov yj] proposed a relatively simple and effective method for generating 
Dicke states in a string of ions that utilizes frequency-chirped optical pulses. The method described 
in Ref. yJ] can be used to generate Dicke states for arbitrary numbers of particles and excitations. 
For small numbers of excitations m, resonant pulses can also be used to produce Dicke states 
although in this case perfect fidelity is possible only when m = 1 . Hume et al. [fioll analyzed the 
fidelity attainable by a method using resonant pulses and experimentally demonstrated it using two 
25]y[g+ ions and an ancillary ^^Al"^ ion. 

In this article we report the generation of an entangled state of two ions, \D^2^) = (1/ V2)(|iT) + 
ITi)), by using rapid adiabatic passage (RAP) [|lll - ll3n in a sideband transition of S 1/2-D5/2 of 
"^^Ca"^. We also analyze the effects of using optical pulses with time-dependent envelopes for RAR 
Such pulses produce time-dependent AC Stark shifts that may result in violation of adiabaticity. 
One advantage of using adiabatic passage to generate Dicke states is its robustness against param- 
eter variations. Using an AC-Stark-shift compensator that cancels time-dependent AC Stark shifts, 
the robustness of entanglement generation is verified by varying the width and peak intensity of 
the optical pulse used for RAR The method employed here can be extended to generate larger 
Dicke states with more particles and excitations. 
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II. PRINCIPLES 



Fig. \T\ depicts the basic concept for generating Dicke states using RAP. The internal states of 
two ions and the motional Fock states of their axial center-of-mass (COM) mode are respectively 
denoted as \s1S2) and \n), where Si(s2) =i, T represents a qubit state of each ion and n denotes the 
motional quantum number. 

The ions are initially prepared in a motional Fock state |1) using either n pulses or RAP 
pulses on the blue (red) sideband and a carrier transition that are addressed to only one of the ions; 
e.g., a blue sideband n pulse followed by a carrier n pulse, both of which are addressed to the first 
ion, to transfer the ions as |0) ITi) |1) ^ lii) |1)- 

An entangling operation is performed with a laser beam tuned closely to ojl = (Oq - ojy, which 
excites the red sideband transition, where ojl (coq) is the laser (atomic resonance) frequency and 
ojv is the secular frequency of the axial COM mode. This preserves the total number of excitations 
Ne = Na + n, where A^^ is the number of atoms in |T); hence, |1) is coupled only to |iT) |0) and 
ITi) |0), and the relevant states form a V-shape diagram. 

When the ions are excited by a single laser, the time dependences of the Rabi frequency for 
each ion are the same up to a proportionality factor. In addition, the detuning for each ion will also 
be the same when the difference in level shifts due to Zeeman or Stark effects can be neglected. In 
such a case, the system can be further simplified by a unitary transformation in the upper manifold 
with n = spanned by |iT) |0) and |Ti) |0), and it becomes equivalent to a two-state system and an 
uncoupled (dark) state (right diagram in Fig. 1(a)). This is one of the simplest examples of a more 



1511 . If the two ions have the 



general prescription known as the Morris-Shore transformation [|14l 
Rabi frequencies of the same time dependence, the two-state system will consist of |1) and 
a Dicke state (or a Bell state) \D^^^) |0) = (1/ V2)(|iT) + IT-i)) |0), with the uncoupled state being 
equal to another Bell state with the same A^,„ (1/ V2)(|iT) - ITi)) |0). 

Thus, starting from |1), either a n pulse or RAP on the red sideband transition can be used 
to transfer the population to the desired final state ^) |0). Fig. 1(b) shows the adiabatic potentials 
for the case of RAP and plots the energy levels with different detunings for the red sideband pulse. 
The population follows the lower curve to be adiabatically transferred from |1) to \D''^-') |0). 
(Here, the Rabi frequencies are assumed to be constant for simplicity, while it is usually varied in 
actual experiments. ) 

This scheme can be extended to more ions N > 2 and more excitations A^^g > 1 to produce 
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I^w''^) I0)- The latter extension requires extending the Morris-Shore transformation to multilevel 
ladders UM; in this case, the initial state |i ... i) \Ne) is coupled to the final state by a ladder with 
Ne + I energy levels. With unitary operations induced by resonant optical pulses, the dynamics 

' 16D, 



in this multilevel diagram never produce a unit population in the final state \D''!!''') |0) [|10L 



whereas when using RAP, such an objective can be accomplished, at least under ideal conditions 
lit. 



III. EXPERIMENTAL PROCEDURE 



The trapping system and laser setup used in the present study are similar to those de- 
scribed in IitI . In brief, we use two "^''Ca^ in a linear trap with center-of-mass frequencies 
of {a)x,a)y,a)J/2n = (2.4,2.2,0.7) MHz. As the qubit states |i) = |5 1/2, = -1/2) and 
\1) = 1^5/2, my = -3/2) are used, and the transition between them is excited by a titanium- 
sapphire laser beam with a wavelength of 729 nm. The envelopes and frequencies of the optical 
pulses used for RAP are controlled by an acousto-optic modulator (AOM) in a double-pass con- 
figuration. Linear frequency chirping of the rf for the AOM is performed using a direct digital 
synthesizer. The pulse envelope is varied by using a mixer with a digital-to-analog converter as 
an input. Both the direct digital synthesizer and the digital-to-analog converter are controlled by a 
field-programmable gate array. 

Individual addressing is required to prepare motional Fock states in advance to the entangling 
operation, which is here done using AC Stark shifts generated by a 854 nm beam that is off" reso- 
nance with the ^D5/2-^P3/2 transition. This beam, which has a l/e^ intensity radius of 15 /im, is 
applied to the two ions at 45 degrees to their axis, with a small offset from their center to produce 
diff'erent shifts for the ions. Typical shifts for the two ions are ~(10, 100) kHz, respectively. 

After Doppler cooling and ground-state cooling of axial COM and stretch modes, average 
quantum numbers (hcom, ^st) ~ (0.06, 0.06) are obtained, respectively. A motional Fock state 
\n = 1) is prepared by applying a blue sideband and a carrier n pulse in series to one ion using the 
individual addressing technique described above. An entangling operation is then performed by 
applying a chirped pulse to the red sideband transitions of both ions. A chirp width of 200 kHz 
typically used; i.e., the detuning is swept from -100 to -I- 100 kHz. 

We also use an AC-Stark-shift compensator in experiments to verify the robustness of Dicke 
state generation (see ^VTl) . To realize this, we use another AOM in a single-pass configuration. 
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Two RF signals with independently variable frequencies and powers, one of which is for direct 
excitation and the other for AC-Stark-shift compensation, are combined and fed to this AOM. The 
output of this AOM is directed into a single-mode optical fiber and then into the above-mentioned 
double-pass AOM system, and it is applied to ions. 

IV. RESULT OF DICKE STATE GENERATION 

The fidelity for generating a Dicke state F = (d^^^|p = (p^^^if + Pn,n)/^ + 

Re(PiT,Ti)' analyzed by performing global rotations and projective measurements. The sum 
of the diagonal elements, pj,T,iT + PTi,Ti estimated by performing projective measurements im- 
mediately after the RAP pulse. The result is shown in Fig. [2ta). The histogram of fluorescence 
counts obtained by using a 397-nm laser beam is plotted. From this, pj,t,j,T + Pu,U - '^•^^ - ^-^^ 
obtained. 

The term containing oflJ'-diagonal elements, Re(p|^;[^|) is estimated by parity measurements after 
applying a n/2 pulse with a variable phase. The parity after applying a n/2 pulse with phase (p is 
expressed as follows. 

n(0) = tr(UR,/2(cpypRn/2((/>)) 

= 2 (Re(pit,Ti) - Re(Pii,TT) cos 2(p + Im(p^,^^) sin 2(p) , 

where ft = |U)(UI + ITTXTTI - liTXiTI - ITiXTil is the parity operator. K,2((f>) represents 
the rotation operation with a n/2 pulse with phase 0. In cases with perfect fidelity, n(0) is unity 
regardless of 0. In realistic cases, the ofl'sets of sinusoidal fits to the measured values give the 
values of 2Re(pi^^fi). 

Fig. [2b) shows parity measurement results. Open circles represent measured values of 11(0) 
and the solid curve is a sinusoidal fit with period n to the measured values. From this, 2Re(p||j|) = 
0.58 + 0.02 is obtained. Combined with the above result for the diagonal matrix elements, the 
fidelity for Dicke state generation is estimated to be F = 0.66 + 0.03. 



can be transferred to (1/ y/2)(\ii) + |TT)) with a n/2 pulse, and parity measurement after 
such a transfer may also give an indication of entanglement. The filled circles in Fig.[2tb) represent 
parity after Dicke state generation followed by two n/2 pulses with the phase of the second varied. 
The sinusoidal oscillation with period n indicates the presence of (1 / V2)(|ii) -I- |TT)) after the first 
n/2 pulse. By combining this result with that for one n/2 pulse described above, we estimate the 
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fidelity for generation of (1/ V2)(|ii) + |TT)) to be 0.64 + 0.04. 



V. EFFECTS OF AC STARK SHIFTS 

The method used here for generating Dicke states using RAP on sideband transitions requires 
pulses with time-dependent envelopes to perform efficient population transfer. As is well known, 
when intense optical pulses are applied to sideband transitions th ey p roduce large AC Stark shifts 
and can cause undesirable eff'ects such as qubit-phase rotations [18]. In the sideband RAP case 
treated here, relevant adiabatic potentials may be deformed in unexpected ways and adiabaticity 
may be violated. This was avoided in the result presented above by selecting sets of parameter val- 
ues such that the dynamics is not greatly affected by diabatic transitions. However, if we consider 
robustness against variation of parameters such as the pulse width and the peak Rabi frequency, 
which is one of the biggest advantages of the adiabatic method described here, we should manage 
AC Stark shifts appropriately. 

We calculated adiabatic potentials and considered adiabaticity in the case of Dicke-state gener- 
ation using sideband RAP. The calculation is based on the following time-dependent Hamiltonian 
in a rotating frame. 

2 

H(t) = ^ IT)y (Tly + ^^^vfl^fl 
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;=i L ^ ^ -I 

Here, l*)^- (5 =i, T) are state vectors for the jth ion, (t+j = |T)^- and &,j = |T)y <il^- + <TI;- 
Q.{t) is the time-dependent Rabi frequency for carrier transitions, and 6{t) is the detuning from 
their resonance. 77 and a are respectively the Lamb-Dicke parameter and the annihilation operator 
for the axial COM mode. 

We chose the following five basis states for the calculation, |U)|0), lii)|l), |£>) |0), |I>)|1), 
and |TT)|0), where \D) represents \D[^^) = (1/ V2)(|iT) + ITi» [see Fig. Ha)]. These consist of 
the two states |1) and \D) |0), which are connected to each other by red sideband RAP, and 
states that are directly connected to these two states by carrier excitation. The laser is assumed 
to illuminate both ions equally. In that case, the antisymmetric Bell state (1/ V2)(|iT) - ILL)) 



becomes uncoupled with other internal states and can be excluded from consideration [|l4l 



The existence of |TT) |1) can be neglected since it is only indirectly coupled with |1) via |Z)) |1) 



through off-resonant carrier couplings. In fact, even |0) and |tT) |0) can be neglected when 
calculating energy shifts, since the AC Stark shifts of \D) |0) due to the couplings with those states 
exactly cancel each other and have negligible effects. Here, we consider these states for the sake 
of completeness. Therefore, the only non-negligible shift originates from the coupling between 
IU)|l)and ID) 11). 

Fig.Hb) shows the time dependence of the Rabi frequency (solid curve) and detuning (dashed 
curve) of the assumed optical pulse. The Rabi frequencies represent those for the carrier transi- 
tion. Fig. Oc) shows the adiabatic potentials for relevant dressed states (solid curves), obtained 
by diagonalizing a Hamiltonian derived from the one given above. The two potentials roughly 
correspond to the bare states |1) and \D) |0), where \D) = {D^ ^Y (The quantum number of the 
optical field is omitted for simplicity.) 

The bare-state potentials for |1) and \D) |0) intersect at the center (about 0.3 ms), where the 
adiabatic potentials for the dressed states approach each other, as expected. However, they also 
approach each other at ~0.1 1 ms, which is not expected in the case of RAP in two-level systems. 
Under some conditions, diabatic transitions could occur at such points, reducing the fidelity. 

The probability for diabatic transitions from state \i(t)) to \j{t)) is estimated to be jiol 



Pi^j < max 



iOjiit) 

where ay;(0 = {j(t)\(d/dt) \i(t)) and ojji(t) = [Ey(0 - £,(0]/^, with 8,(0 being the adiabatic potential 
for \i(t)). In short, the probability is determined by the maximum value of \aji{t)/a)ji{t)\^ . Fig. 
[3td) shows the calculated values of \a ji(t) / co ji(t)\^ . At points near 0.11 ms, where the adiabatic 
potentials approach each other, this value has a local maximum, which indicates that adiabaticity 
is about to be violated. At slightly higher intensities, the probability of diabatic transitions may 
become fatally large. 

For comparison, potentials for an ideal case when the carrier couplings are omitted are also 
calculated. The dashed curves in Fig. ^c) are calculated potentials for when the matrix elements 
for all the carrier couplings (|ii) |0)-|D) |0), |TT) \0)-\D) |0), |ii) |1)-|D) |1)) are set to be zero. 
These potentials correspond to the case in which the energy shifts due to the carrier couplings are 
perfectly compensated by the use of an AC-Stark-shift compensator (see the next section). They 
have clean shapes and there is only one point where they approach each other. Correspondingly, 
as shown by the dashed curve in Fig. Od), \a pit) / oj ji{t)\^ peaks only at the center, where the 
adiabaticity is easily satisfied since the Rabi frequency is a maximum at that point. 
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VI. VERIFICATION OF ROBUSTNESS OF DICKE STATE GENERATION 



Adiabaticity can be maintained by performing AC-Stark-shift compensation [llSll . This can be 
done by applying another optical pulse such that the AC Stark shift that it produces nearly cancels 
the original shift. Using the AC-Stark-shift compensator described in ^Illl we experimentally 
verified the robustness of Dicke state generation against parameter variations. The compensator is 
used during both the blue sideband pulse in the Fock state preparation process and the red sideband 
RAP pulse. It is detuned in the opposite direction to the main beam for direct excitation in each 
period, and the absolute value of its detuning is set around a value determined by the power ratio 
of the two beams such that the two shifts nearly cancel each other. Here, the optical power of the 
compensator is set to ~60% of the main beam; accordingly, the absolute value of the detuning is 
centered on ~ 400 kHz. 

The fidelity measurement results for when the peak Rabi frequency and the pulse width are 
varied are indicated by the circles in Fig. St a) and (b), respectively. In both cases, the fidelity 
exceeds 0.5 over wide variations in the parameters (0.5 represents the inseparability criterion for 
two-particle entangled states 1I20I1 ). 



VII. DISCUSSION OF FIDELITY 



There are several possible causes for the fidelity of Dicke state generation being limited to 
below 0.66. The error consists mainly of that produced when preparing motional Fock states and 
in the entangling operation with RAR The causes for the former include insufficient cooling to the 
motional ground state and an addressing error, and the causes of the latter include phase relaxation 
processes due to fluctuations in the laser frequency and magnetic field and fluctuations in residual 
AC Stark shifts due to beam jitter. 

To analyze such infidelity factors, the values of pj,t,iT +PTi,Ti '^^^(Pitti)^ which respectively 
reflect the diagonal and off-diagonal density matrix elements, are shown in Fig. H] (as crosses and 
asterisks, respectively). Relatively high values are obtained for the former, whereas the latter 
is limited to below 0.5 over both parameter ranges and it decays with increasing pulse width. 
Based on these facts, we speculate that the error generated in preparing motional Fock states is not 
significant, since otherwise we would observe a large reduction in the diagonal matrix elements. 
We can therefore conclude that the infidelity factors are mainly associated with the entangling 
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operation using RAP. We speculate that phase relaxation processes during the entangling operation 
mainly cause the infidelity, which is consistent with the decay with increasing pulse width. We 
expect that the fidelity can be improved by using resonant pulses with improved 729-nm laser 
linewidths (currently <400 Hz) and using transitions that are less sensitive to the magnetic field, 
such as 5i/2(mj = -Il2)-Dsi2{mj, = -1/2). 

VIII. CONCLUSION 

We have demonstrated generation of entangled states using an adiabatic method and demon- 
strated the robustness of the generation process against parameter variations. The method de- 
scribed here can be extended to larger numbers of particles and excitations. Recently, Wieczorek 
et al. [6] have shown experimentally that entangled states of inequivalent classes are obtained from 
six-photon Dicke states l^^g ^) via projection measurements on a few of their qubits. In a similar 
way, large atomic Dicke states generated by methods similar to that described here may be used 
as entanglement resources for quantum information processing applications. 
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FIG. 1: (Color online) Basic concept for generating Dicke states using rapid adiabatic passage, (a) (left) 
Relevant states. A motional Fock state |1) is coupled to |J,t) |0) and |tJ,) |0) with a red sideband pulse, 
(right) After applying a basis transformation (Morris-Shore transformation) to the upper two states, the 
system is decomposed into a two-state system and an uncoupled (dark) state, (b) Dressed state potentials 
for the two-state system plotted against the detuning of the laser from the red sideband resonance. In 
rapid adiabatic passage, the ions move from the left along the potential (indicated by the solid curve) to 
1^2 ''} |0) = (1/ V2)(|J,t) -I- ITJ-)) |0). Here, the Rabi frequencies are assumed to be constant for simplicity, 



while it is usually varied in actual experiments. 
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FIG. 2: (Color online) (a) Result of projection (population) measurement after Dicke state generation. A 
histogram of fluorescence counts obtained by shining a 397-nm laser beam for 7 ms after the generation 
of a Dicke state is shown. The left, middle, and right peaks (around 0, 70, and 140 counts, respectively) 
correspond to Pii,Pii+Pii, and P^, respectively, where P^, 52(^1, ^2 =i, t) represents the population in the 
respective basis state. From this, pj,-f,j,-f + Pti,Ti - 0-74 + 0.06 is obtained, (b) Result of parity measurement 
after Dicke state generation. The measured parity after Dicke state generation followed by the application of 
a n/2 pulse is plotted against the phase of the n/2 pulse phase (open circles). The solid curve is a sinusoidal 
fit with period k to the measured values. From this, 2Re(pj^-[- -[-i) = 0.58 ± 0.02 is obtained. Filled circles 
represent parities after Dicke state generation followed by two n/2 pulses and the phase of the second pulse 
is variable. The sinusoidal oscillation with period jt indicates the presence of (1/ V2)(|J,J,) + |tT)) after the 
first n/2 pulse. 
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FIG. 3: (Color online) Effects of AC Stark shifts due to optical pulses for rapid adiabatic passage, (a) 
Energy levels used for the calculation (horizontal solid lines). |D) represents ^) = (1/ V2)(|J,t) + ITi))- 
The bold arrow indicates the transition used for rapid adiabatic passage and the other arrows represent 
couplings considered here for calculating the AC Stark shifts, (b) Variation of the envelope (solid curve) 
and detuning (dashed line) of an optical pulse for rapid adiabatic passage against time. The envelope is 
assumed to be a Gaussian function, and the detuning is assumed to be linear, (c) Variation of adiabatic 
energies of relevant levels against time (solid curves). They approach each other near 0.11 ms, as well 
as near 0.3 ms. The dashed curves are shown for comparison; they represent adiabatic energies for an 
ideal case in which Hamiltonian matrix elements for carrier transitions ( |1)-|D) |1), |D) |0)-|J.J,) |0), 
I'D) |0>-|tT) |0» are set to be zero. These potentials only approach each other near 0.3 ms. (d) |ay;(0/^^yr(Of 
indicating violation of adiabaticity is plotted (see text). 
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FIG. 4: (Color online) Verification of robustness in generation of Dicke states using rapid adiabatic passage. 
Gaussian pulses producing time-dependent Rabi frequency 0(0 = Opeak exp[-f^/2cr^] are applied as RAP 
pulses and values of obtained fidelity (circles) are estimated, with (a) the pulse width Icr and (b) the peak 
Rabi frequency Qpeak;/27r varied. Here, the Rabi frequency is for the carrier transition. Qpeak/27r is fixed to 
145 kHz in (a) and 2cr is fixed to 244 //s in (b). The total duration of each RAP pulse is 2.36 times the pulse 
width 2o", and during the pulse duration the excitation laser is swept from -100 kHz below to 100 kHz 
above the resonance of the red sideband transition. The error bars represents 68% confidence intervals. The 
horizontal dashed lines at a fidelity of 0.5 represent the threshold for inseparable states (see text). Crosses 
and asterisks indicate values of pj,-|-,j,T + PTi.Ti ^^'^ 2Re(p|-|- ^j^), respectively, which are used for discussions 
about infidelity in gVIII 
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